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ABSTRACT 

We present the detection and characterization of a peculiar low-mass protostar (IRAS 22129+7000) 
located ~0.4pc from Cederblad 201 Photodissociation Region (PDR) and ~0.2pc from the HH450 jet. 
The cold circumstellar envelope surrounding the object has been mapped through its 1.2 mm dust 
continuum emission with IRAM-SOm /MAMBO. The deeply embedded protostar is clearly detected 
with Spitzer/MlPS (70 (jem), IRS (20-35 fim) and IRAC (4.5, 5.8, and 8 fan) but also in the K s band 
(2.15 /xm). Given the large near- and mid-IR excess in its spectral energy distribution, but large 
submillimeter-to-bolometric luminosity ratio (~2%), IRAS 22129+7000 must be a transition Class 0/1 
source and/or a multiple stellar system. Targeted observations of several molecular lines from CO, 
13 CO, C 18 0, HCO + and DCO + have been obtained. The presence of a collimated molecular outflow 
mapped with the CSO telescope in the CO J=3-2 line suggests that the protostar/disk system is still 
accreting material from its natal envelope. Indeed, optically thick line profiles from high density tracers 
such as HCO + J=l— show a red-shifted-absorption asymmetry reminiscent of inward motions. We 
construct a preliminary physical model of the circumstellar envelope (including radial density and 
temperature gradients, velocity field and turbulence) that reproduces the observed line profiles and 
estimates the ionization fraction. The presence of both mechanical and (non-ionizing) FUV-radiative 
input makes the region an interesting case to study triggered star formation. 

Subject headings: ISM : individual (IRAS 22129+7000; Ced201; B175/L1219) — ISM : jets and 
outflows — stars: formation — stars: pre-main sequence 



1. INTRODUCTION 

Low-mass stars (sunlike; M*<8Afo) form from the 
gravitational collapse of dense molecular cores within gi- 
ant molecular clouds. During the earliest stages of evolu- 
tion, embryonic protostars are deeply embedded in cold 
and dusty envelopes of infalling material. In recent years, 
large millimeter (mm) and submm continuum surveys 
have been carried out to detect the cold dust emission 
associated with these early birth-sites of star formation 
(Andre et al. 1993; Johnstone et al. 2000). During col- 
lapse, conservation of angular momentum combined with 
infall along the magnetic field lines leads to the forma- 
tion of a rotating protoplanetary disk that drives the 
accretion process. At the same time, both mass and an- 
gular momentum are removed from the system by the 
onset of jets, collimated flows and the magnetic braking 
action (Bontemps et al. 1996; Cabrit et al. 2007). The 
resulting molecular outflow starts to erode and sweep up 
part of the natal envelope, contributing for the clear- 
ing of the circumstellar material and the termination of 
the infall phase (Arce & Sargent 2006). The so-called 
Class sources are deeply embedded protostars at their 

1 Based on observations obtained with the CSO telescope. The 
CSO telescope is operated by the California Institute of Technology 
under funding from the NSF, Grant No. AST-0540882. 

2 Observations also obtained with IRAM telescopes, supported 
by INSU/CNRS (France), MPG (Germany), and IGN (Spain). 

3 This work is based on observations made with the Spitzcr Space 
Telescope, which is operated by the Jet Propulsion Laboratory, 
California Institute of Technology under a contract with NASA. 



main accretion phase and represent the earliest type 
of young stellar objects (YSOs). Observationally they 
should show: (i) Evidence of a central YSO (as revealed 
by the presence of a collimated outflow or an internal 
heating source); (ii) Centrally peaked but extended mas- 
sive envelope (as traced by mm and submm dust contin- 
uum emission); (iii) High submm (A>350/zm) to bolo- 
metric luminosity ratio (L S7nm /L(, ; >0.5%), and a spec- 
tral energy distribution (SED) similar to a single gray- 
body at 15-30 K (Andre et al. 2000). On the other hand, 
Class I sources correspond to the late accretion phase 
(roughly half of the initial envelope mass is accreted in 
the Class stage), they don't satisfy (iii) and can easily 
be detected at near-IR wavelengths (Lada 1999). 

The detection of YSOs that do not exactly match the 
above canonical categories can provide the missing pieces 
of our current understanding of early disk evolution and 
early dissipation of the circumstellar envelope. These 
sources can thus potentially help to address, for example, 
how and when the velocity field around YSOs changes 
from infall to rotation dominated, which fraction of the 
initial mass envelope is dissipated and not accreted onto 
the star/disk, as well as the timescales for such dissipa- 
tion (Brinch et al. 2007; J0rgensen et al. 2007) 

In this paper we present the detection of a low- 
mass protostar in the Bok globule B175 (L1219) within 
the Cepheus Flare, a nearby (^400 pc) star form- 
ing region ~85pc above the galactic plane (Kun 1998; 
Nikolic & Kun 2004). This cloud hosts the enig- 
matic Cederblad 201 reflection nebula (Wolf 1908; 
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TABLE 1 
-7000 Photometric data 



v LSR [km s ] v LSR [km s ] 

Fig. 1. — Middle panel: IRAM-30m/MAMBO 1.2 mm-continuum emission of the region. The image is centered at IRAS 22129+7000. 
White contours are displayed from 22 (2cr) to 55 mjy/15"-beam in steps of 11 mjy/15"-beam (lcr). The CO J=3-2 integrated intensity 
(from 5 to 9 K km s _1 in steps of 1 K km s — 1 ) from the outflow is shown in red (from v=-4 to -2 km s — 1 ) and blue (from v=-6 to -8 km s -1 ). 
The integrated CO J=3-2 line intensity from Ced 201 PDR (v=-5.5 to -4.5 km s - 1 ) is shown in black contours (from 6 to 20 K km s -1 in 
steps of 1 K km s" 1 ). Left panels: 12 CO J=2-l and 4-3, HCO+ J=l-0 and DCO+ J=3-2 lines towards IRAS 22129+7000 position. Right 
panels: Declination raster in the HCO + J=l-0 line. The vertical dashed line marks the rest velocity of the cloud (— 5kms _1 ). 

Cesarsky et al. 2000), probably the result of a chance 
encounter with a B9.5V runaway star (BD+69°1231, 
T e ff ~10,000 K), which is moving through the cloud at 
~12kms~ 1 (Witt et al. 1987). The star is thus not a 
only a source of dissociating photons (the FUV radiation 
field is ~200 times the mean interstellar radiation field) 
but also a source of shocks and turbulence in the region. 
This interaction enables us to probe matter in a different 
situation than in classical PDRs in which the distance 
from the exciting star to the illuminated regions is large 
and/or in which the star was formed locally. In addition, 
the extended cloud hosts a Herbig-Haro jet (HH450), 
and several arsec-scale filaments of the supernova rem- 
nant (G110.3+11.3) that are rapidly approaching to the 
region (Bally & Reipurth 2001). The presence of both 
mechanical (from shocks and outflows) and radiative in- 
put (from FUV photons) makes the region an interesting 
dynamics laboratory to study how star formation is pro- 
ceeding through the cloud. 

2. OBSERVATIONS 

A 1.2 mm continuum emission map of the region was 
obtained at the IRAM-30m (Pico Veleta, Spain) in 
March 2005 using the 117-channel MAMBO bolome- 
ter array. The angular resolution is ^11". A 
fast-mapping mode was used to map the region 
(Teyssier & Sievers 1999). The total integration time 
was lhour, achieving a rms noise of 7mJy/ll"-beam. 
Sky noise substraction and data analysis were carried out 
with the M0PSIC software (Zylka 1998). A 31"x22" full- 
width at half-maximum (FWHM) condensation with a 
peak intensity of 60±7mJy/ll"-beam is detected ^225" 
East of Ced 201 PDR. The integrated flux is 249±50mJy, 
which reflects the estimated absolute calibration uncer- 
tainty. The resulting 1.2 mm emission map is shown in 
Figure 1 (smoothed to a beam resolution of 15"). 

The HCO + J=l-0 line was observed also at the -30m 
telescope during September 2005 using the B100 single 
sideband receiver with a channel resolution of 40 kHz. 
The main beam efficiency at 89 GHz is 0.82 and the 



A 


Flux 


Telescope 


Aperture 


(H 


(Jy) 






1.25 


t3.20xl0 -4 


2MASS 


3" 


1.65 


t3.60xl0" 4 


2MASS 


3" 


2.15 


3.00X10" 3 


2MASS 


3" 


4.5 


0.17 


Spitzer/IRAC 


12.2" 


5.8 


0.38 


Spitzer/IRAC 


12.2" 


8.0 


0.31 


Spitzer/IRAC 


12.2" 


12 


+0.39 


IRAS 


5' 


25 


1.13 


IRAS 


5' 


60 


3.46 


IRAS 


5' 


70 


6.10 


Spitzer/MIPS 


1' 


100 


16.00 


IRAS 


5' 


200 


22.00 


KAO 


~1' 


1200 


0.25 


IRAM-30m/MAMBO 


31"x22"* 



Note. 



' Upper limit. 'Source FWHM size given. 



angular resolution is 28". The line was observed to- 
wards the dust condensation and in a declination raster 
crossing Ced 201 PDR (see Fig. 1). Additional observa- 
tions were performed at the CSO telescope (Mauna Kea, 
Hawaii) during July 2006 and June 2007. The telescope 
is equipped with SIS receivers operated in double side- 
band mode. The CO J=3-2 map (Figs. 1 and 2) was 
made "on-the-fly" . The DCO + J=3-2 line was mapped 
in a ~40" x30" region around the dust condensation, 
whereas pointed observations in the CO J=2-l and 4-3 
lines were taken towards the dust emission peak (Fig. 1). 
The spectra were analyzed with a 1024 channel acousto- 
optic spectrometer with a total bandwidth of 50 MHz and 
a resolution of ~100kHz. The main beam efficiencies of 
the telescope were 0.70, 0.75 and 0.53 at 230, 345, and 
461 GHz respectively. The angular resolution is ~30" at 



CO J=2-l, ~20" at CO J=3-2 and ■ 
Data was processed with GILDAS. 



45" at CO J=4-3. 
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Fig. 2. — Spitzer /IRAC 8 fim image of Ced201 region including the detection of IRAS 22129+7000 protostar. The bright reflexion nebula 
(the PDR) is clearly seen at the west of the region. A well collimated outflow (as traced by CO J=3— 2 high velocity line wing emission) 
emerges from the protostar position. Contours show the integrated intensity (from 5 to 9 K km s _1 in steps of 1 K km s" 1 ) in the red 
(from w=-4 to -2 km s _1 ) and blue (from v=-6 to -8 km s _1 ) line wings respectively. 



The Ced 201 region was also observed with Spitzer 
as part of our SPECPDR program 4 (Joblin et al. 2005) 
with the Infrared Array Camera (IRAC) at 4.5, 5.8, and 
8 /jm and with the Infrared Spectrograph (IRS). IRS 
data were obtained in the spectral mapping mode with 
the LL1 module (spectral range: 20 to 35 /jm). The spec- 
tral cube of the source was assembled using the CUBISM 
software (Smith et al. 2007) from the Basic Calibrated 
Data (BCD) files. The 20 to 35 /xm data was obtained 
by integrating the cube spatially (see Berne et al. 2007) . 
The 70 /xm Multiband Imaging Photometer for Spitzer 
(MIPS) data that are part of the PID 30790 program 
(Fazio et al.) were retrieved from the Spitzer archive. 
MIPS and IRAC post-BCD calibrated images were used 
directly without additional processing. 

3. RESULTS: DETECTION OF A YOUNG PROTOSTAR 

Figure 2 shows the detection of a compact contin- 
uum source at ~8/zm (warm dust) located ~0.4pc East 
of Ced 201 PDR (Berne et al. 2007). The compact 
object is also detected at 4.5, 5.8, 20-35 and 70 ^m 
with Spitzer. The near-IR source position agrees with 
the continuum peak of the 1.2 mm extended emission 
(cold dust) detected in a larger field-of-view (Fig. 1). 
The best coordinates inferred from IRAC images are 
a 200 o = 22 h 14 m 08.3 s , <5 2 ooo = +70°15'06.7", which coin- 
cides with the IRAS 22129+7000 source, clearly a YSO. 
A well collimated molecular outflow emanating from the 
object is seen along a SE-NW axis in the CO J=3-2 line 
(Figs. 2 and 1). This outflow was previously detected 
in CO J=l-0 at lower angular resolution by Nikolic & 
Kun (2004), who first suggested that IRAS 22129+7000 
could be the driving source. We now detect the near- 
and mid-IR scattered light from the embedded proto- 
star and the mm thermal emission from the surrounding 
dusty cocoon. In the following we examine the nature of 
this source and whether or not the presence of a PDR 
plays a role in the star formation process in the region. 

3.1. Peculiar SED Characteristics 

In order to investigate the emitted power in the differ- 
ent wavelength domains, the SED of IRAS 22129+7000 
has been built with our MAMBO and Spitzer ob- 
servations and complemented with previous 2MASS 

4 http://www.cesr.fr/~joblin/SPECPDR_public/Home 



(2.15/mi), IRAS (12, 25, 60 and 100 fim) and KAO 
(200 ^m; Casey 1991) detections (Table 1). We fit- 
ted the SED with two graybodies, B\(T) (1 - e" TA )fi, 
where f2 is the solid angle subtended by the emit- 
ting region. In the far-IR and mm domain t\ oc 
K1200 (1200/A) and T can be associated with the dust 
opacity and color temperature (— T^) of the extended 
envelope. Here we take a dust opacity (per gas+dust 
mass column density) of ki 2 oo=0.01 cm 2 g _1 at 1200 /mi, 
the usual value adopted for the dust in protostellar 
envelopes (Andre et al. 2000). Because of the implic- 
itly assumed grain growth and ice mantle formation 
(Ossenkopf & Henning 1994), this value is a factor ^5 
larger than K1200 for standard interstellar grains. Q is 
fixed as the FWHM size of the envelope inferred from the 
1.2 mm continuum map, while the column density of ma- 
terial is varied between Ay ~5 and 30. Satisfactory fits 
are obtained for Td=24-19 K and dust spectral indexes 
^=1.4-1.8 respectively (see the best fit as continuous blue 
curve in Fig. 3). A luminosity of Lb o /=5.5±O.5L and a 
"bolometric temperature" of Tboi — 175K (the tempera- 
ture of a blackbody with the same SED mean frequency; 
Myers & Ladd 1993) are inferred from the full fit. 

In a first inspection, the resulting SED (Fig. 3) does 
not look like a typical Class source in the near- and 
mid-IR (see e.g., Whitney et al. 2003). The improved 
sensitivity of the IR instrumentation compared to that 
available in the 90's has allowed the detection of scat- 
tered light from several Class sources below ~T0 /im 
(Tobin et al. 2007). In addition, the K s -band emission 
due to scattered light observed through a cavity opened 
by the outflow (with a favorable inclination angle) and/or 
H 2 rovibrational emission from the outflow base, has also 
been detected in a few Class sources (e.g., Tachihara 
et al. 2007). IRAS 22129+7000 shows these observable 
characteristics but it also displays an intriguing strong, 
and rather flat, near- and mid-IR emission (as shown by 
the relatively high Th /). In fact, Spitzer fluxes towards 
IRAS 22129+7000 are consistent with the typical val- 
ues observed in Class I sources with the same luminosity 
(Furlan et al. 2008), but an order of magnitude higher 
than in Class sources (Tobin et al. 2007). On the other 
hand, IRAS 22129+7000 shows a very large fraction of 
submm (A >350 /im) to bolometric luminosity (~2%; a 
signature of younger Class sources still surrounded by 
massive dusty envelopes). 
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Fig. 3.— SED of IRAS 22129+7000 combining data in the H 
and K s bands (2MASS; Skrutskie et al. 2006); 4.5, 5.8 and 8/fln 
(Spitzer/IRAC); 20-35 /jm (Spitzer/IRS, smaller points), 12, 25, 60 
and 100 fun (IRAS), 70 /an (Spitzer/MIPS), 200 Aim (KAO: Casey 
1991) and 1200 jim (MAMBO). The Spitzer/IRS fluxes towards 
the source agree with those measured by IRAS at 25/xm probing 
the compact nature of the mid— IR emission. The dashed curve is 
a simple fit using two graybodies with different temperatures. The 
continuous blue curve represents the circumstellar envelope with 
T [J =22K and a spectral index of 1.6. The magenta pointed— and 
pointed-dashed curves show the best 2D radiative transfer models 
obtained with the Robitaille et al. (2007) fitter for two different 
inclinations (i = corresponds to "pole-on" and i = 90 to "edge- 
on" configurations). 

In an attempt to better understand the peculiar IR 
shape we tried to fit the SED with the 2D radiative 
transfer tool by Robitaille et al. (2007). The best mod- 
els fitting the data are shown in Fig. 3 for two different 
inclination angles. While the long-wavelength emission 
from the circumstellar envelope is well reproduced by the 
models (and is almost not affected by the viewing angle), 
the short- wavelength emission from the protostellar/disk 
system is clearly underestimated and strongly depends 
on inclination. Therefore, a small viewing angle favors 
the detection of near- and mid-IR scattered light from 
the embedded protostar through the cavities cleared by 
the outflow. However, the blue- and red-shifted lobes 
of the CO outflow (Fig. 2) appear relatively well sepa- 
rated, and thus the source can not be completely "pole- 
in" (Cabrit & Bertout 1986). In summary, the observed 
SED is peculiar, with a near- and mid-IR excess that 
can not be fully explained in terms of inclination effects. 

3.2. Circumstellar Envelope: Radial Analysis 

The cold dust surrounding the YSO has been mapped 
through its 1.2mm emission (Fig. 1). As mm continuum 
emission is generally optically thin, the measured fluxes 
are proportional to the temperature-weighted mass con- 
tent assuming that dust grains properties are known. In 
order to trace the mass distribution in the envelope, the 
radial intensity profile I(r) has been analyzed by aver- 
aging the 1.2 mm fluxes within equidistant circular an- 
nulus (Fig. 4). An envelope radius of R ou t —15,000 AU 
is inferred. Assuming spherical symmetry, the intensity 
profile I(r) oc r~ m contains information on the density 
distribution, p(r) oc r~ p , and temperature distribution, 
Td(r) oc r~ q , with p — m + 1 — q (e.g., Adams 1991). 

Before convolving the expected intensity profile with 
the bolometer angular resolution, we divide the enve- 
lope in two regions (inner and outer) and try to fit the 



TABLE 2 

IRAS 22129+700: SED PARAMETERS 





Tboi 








(L©) 


(K) 


CM©) 


(%) c 


(M /L Q ) 


5.5+0.5 


~175 


0.2 a -0.7 6 


~2 


0.04 a -0.13 6 



a Assuming T d =20K, ki. 2 =0.01 cm 2 g~ x and r <5,000AU. 
b Mass within r=R out <15,000 AU. 

c Lstti?ti is the "submm" luminosity radiated at A >350 /^m. 



observed 1.2 mm emission profile with two power-laws. 
Fig. 4 shows that the dust emission in the outer envelope 
(r>5,000AU >one beam) follows a very steep profile, 
with m ~1.9, connecting the envelope with the lower den- 
sity ambient cloud at r > R ou t ■ Since the outer envelope 
angular radii are larger than the beam's FWHM, the con- 
volution does not modify the slope of the intensity profile 
(e.g., Adams 1991; Motte & Andre 2001). For the same 
reason, our observations can not resolve accurately the 
slope of the inner envelope (r < R in ~5,000 AU). We take 
m<0.5 as the expected value in a free-falling envelope. 
The line radiative transfer models shown in Sect. 3.4 are 
not inconsistent with this assumption. 

Taking again a dust opacity of Ki2oo=0.01 cm 2 g _1 and 
T^=20K (consistent with the SED analysis) we com- 
pute a total circumstellar mass of M en „=0.7MQ within 
t < Rout- However, it is unlikely that all this mass can 
be finally accreted by the protostar/disk. First because 
the outflow will disperse a large fraction of it, and sec- 
ond, because differential rotation in the envelope is ex- 
pected to decouple the inner rapidly infalling and rotat- 
ing layers from the slowly collapsing and rotating outer 
regions (Belloche et al. 2002). Hence, we also computed 
M™„=O.2M 01 the mass within one 11 "-beam (r<Ri n ) 
as representative of the mass inside the inner envelope. 
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Fig. 4. — 1.2 mm radial intensity profile around 
IRAS 22129+7000. The intensity (in mjy/ll"-beam) has been av- 
eraged over circular annulus in steps of 3.5". The error bars show 
the standard deviation within each annulus. The dashed line repre- 
sent a Gaussian beam with a FWHM of 11". The right axis shows 
an estimate of the H2 column density assuming ki,2=0.01 cm 2 g — 1 
and T ( j=20 K. The red curve shows a fit to the data including 
2 intensity power-law profiles (convolved with the beam), with 
Z(r) oc r" 1 ' 9 for r >5,000AU and I(r) oc r"°' 5 for r <5,000AU. 
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The steep slope of the 1.2 mm intensity profile suggests 
that a temperature gradient is also present. In particu- 
lar, the dust thermal profile through the inner envelope 
should be dominated by accretion heating from the pro- 
tostar/disk system with (Terebey et al. 1993): 

TAr) » 38 k( -—)~ Q ( h. V' M) 

v ; V100AU/ \L Q J 

where q depends on the dust grain properties through 
q = 2/(4 + /?). Given the proximity of Ced 201 PDR, 
FUV photons from BD+69°1231 star also heat the outer 
layers of IRAS 22129+7000 envelope. Penetration of this 
external FUV radiation field into the envelope greatly de- 
pends on the scattering properties of the dust particles, 
with larger penetration depths as grains grow toward big- 
ger grains (Goicoechea & Le Bourlot 2007) . In the follow- 
ing we adopt a temperature profile that follows Eq. (1), 
reaches a minimum of ~10 K at r —10,000 AU and then 
increases up to ~15 K at R ou t due to external FUV heat- 
ing (this temperature is consistent with our analysis of 
CO lines in the ambient cloud). Taking q=0.36 (using 
/3=1.6 from the SED) and assuming L± ~ L^oi, the es- 
timated mass of the envelope increases by a factor —2 
respect to the isothermal case because Td{r) drops be- 
low 20 K in a large fraction of the envelope. This differ- 
ence adds to the dominant source of uncertainty in the 
mass determination, the dust composition and associated 
opacity. 

Taking into account these caveats, the observed 1.2mm 
intensity profile in the outer envelope is equivalent to a 
steep density profile p(r > R in ) oc r~ 2 ' 54 with a number 
density of n(H 2 )~(l-2) x 10 5 cm~ 3 at R in . Under these 
assumptions, the (beam averaged) extinction towards the 
peak (e.g., Motte & Andre 2001) is Ay ~20-30, a factor 
~10 larger than the value implied by our 1.2 mm contin- 
uum observations in the ambient cloud. 

3.3. Gas Kinematics: Outflow and Infall? 

Figure 2 shows the outflow emanating from 
IRAS 22129+7000 and extending -0.3 pc in the NW di- 
rection (redshifted lobe) and ~0.2 pc in the SE direction 
(blueshifted lobe). In addition, IRAS 22129+7000 is the 
driving source of the Herbig Haro HH450 jet detected 
-0.2 pc SE from the object (Bally & Reipurth 2001). 
The CO outflow is well collimated, with an aspect 
ratio of ~3-4 (the projected FWHM flow length/width 
ratio). CO J=2-l to 4-3 line profiles show high velocity 
wing emission from accelerated gas. Typical full line 
widths are consistent with an expansion velocity of 
u oxp ~7kms _1 (Fig. 1). We thus estimate that the dy- 
namical age of the outflow (~ length/w oxp ) is — 35,000yr 
(up to ~10 5 yr if the inclination angle is i=18°). The 
presence of a molecular outflow suggests that the 
protostar/disk system is still accreting material from the 
natal envelope. The high density tracer HCO + J=l-0 
line, with a critical density of nj~°~2 x 10 5 cm -3 , shows 
a deep red-shifted-absorption asymmetry in its optically 
thick line profile (when n(H2)>n^~ J the associated 
J level is predominantly populated by collisions with 
H 2 molecules and T ex , the transition excitation tem- 
perature, tends to the kinetic temperature). The depth 
and velocity position of the absorption dip are due to a 
decrease of T ex at the outer edges of a (likely) infalling 



envelope (left insets in Fig. 1). As low-J CO lines 
easily thermalize at lower densities (n cr ~10 3 ~ 4 cm -3 ), 
they may not be reliable tracers of the presence of 
infall motions but of the varying physical conditions 
in the beam. Among the observed CO lines, only the 
J=4-3 line starts to trace denser (n 4 ~ 3 ~5 x 10 4 cm~ 3 ), 
possibly infalling, gas. In anycase, the HCO + J=l-0 
self-absorption dip is redshifted 0.3 km s -1 from the rest 
velocity of the cloud (— 5kms _1 ) and is not observed 
at nearby positions of the cloud where the line shows a 
pure Gaussian emission profile (right insets in Fig. 1). 
Although higher-J optically thick lines from abundant 
high dipolc moment molecules should be mapped to 
draw definitive conclusions, in the following we as- 
sume that the observed deep HCO + self- absorption is 
produced by inward motions. 

The densest gas in the envelope has been traced with 
the DCO+ J=3-2 line (n 3 r ~ 2 ~2 x 10 6 cnr 3 ) which shows 
compact (<40") emission around IRAS 22129+7000. 
Compared to the broad line wings from the CO out- 
flow, the DCO + J=3-2 line shows the narrowest ob- 
served linewidth (AvFHWM=0.60±0.04kms~ 1 ). On the 
other hand, we used the C 18 J =1-0 line to trace lower 
density gas (e.g., the outer envelope). This line dis- 
plays a broader linewidth (Aufhwm=0. 96+0.02 km s _1 ). 
The radiative transfer calculations presented in the next 
section confirm that both lines are optically thin and 
thus they are not broaden by line saturation. There- 
fore, the different velocity dispersions reflect the distinct 
regions where both line profiles are formed. Assuming 
T fc =10 to 20 K, the observed DCO+ linewidths are 5 
to 3.5 times larger than the expected thermal broad- 
ening. The required non thermal velocity dispersion 
to fit the observed DCO+ J=3-2 linewidth (fi 'om tur- 
bulence, infall or rotation) is a n th — 0.25 km s -1 (with 
A«fhwm = 2.355 x a). 

3.4. Line Radiative Transfer Analysis 

In order to reproduce the observed molecular lines with 
a realistic model of the envelope, we have used a nonlo- 
cal and non-LTE excitation and radiative transfer code 
(Goicoechea et al. 2006). The code models an infalling 
spherical envelope and accounts for line trapping, col- 
lisional excitation 5 , and radiative excitation. Arbitrary 
velocity fields as well as density, temperature, and abun- 
dance radial gradients are included as input parameters. 
This is done by numerical discretization of the envelope 
in spherical shells. The nonlocal radiative transfer prob- 
lem is then simulated by the emission of a determined 
number of model photons (cosmic mm background, con- 
tinuum and line photons) using the Monte Carlo ap- 
proach (Bernes 1979). The steady state statistical equi- 
librium equations are then solved iteratively and the non- 
LTE level populations are determined in the shells. The 
emergent line intensities along each line of sight are fi- 
nally convolved with the telescope angular resolution. 

3.4.1. A First Model for the Circumstellar Envelope 

As our first attempt to understand IRAS 22129+7000 
physical structure, we have modelled a spherical envelope 

CO (and isotopologucs) rotationally inelastic collisional rates 
have been derived from those of Flower (2001). HCO+ and DCO+ 
collisional rates have been derived from those of Flower (1999). 
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Fig. 5.— (a) HCO+ J=l-0, (6) DCO+ J=3-2 and (c) C 18 J=l-0 observations (histograms) towards IRAS 22129+7000. Best radiative 
transfer models for the circumstellar envelope are also shown (red curves). Predicted line profiles have been convolved with the telescope 
angular resolution at each frequency. Intensity scale is in main beam temperature, (d) Density and temperature profiles used in the model. 



with a radius i? oli t=15,000 AU. We take the density and 
temperature profiles (p(r) oc r~ 2 - 54 and T(r) oc r~ - 34 
plus external FUV heating) derived in the previous sec- 
tion for the outer envelope (r > R in ~5,000 AU). Given 
the high densities expected in most of the envelope, we 
assume that gas and dust are thermalized (i.e., Tk=Td). 

Using the above density and temperature structure we 
run the model globally for HCO+, DCO+ and C 18 
and try to find iteratively the best velocity fields and 
abundances fitting the observed line profiles. In or- 
der to reproduce the observed velocity of the CO J=4- 
3 and HCO + J— 1-0 redshifted absorption dips (0.2 
and 0.3 km s -1 ) we adopt a uniform infalling velocity of 
t>inf=0.2kms in the outer envelope. For r < Ri n (the 
inner envelope) we assume p(r) oc r -15 , and v- ul [(r) = 
0.2 (r / Rin)~ 5 i.e., the density and velocity profiles in a 
frccfall inside-out collapsing envelope (Shu et al. 1987). 
Numerical simulations of collapse induced by exter- 
nal compression waves also predict those approximate 
p(r) and v- ln {(r) profiles behind the expanding collapse- 
wavefront (e.g., Hennebelle et al. 2003). These two 
regimes of infall (frccfall and uniform) have been inferred 
in other YSOs with clear evidences of collapse (Belloche 
et al. 2002). Under the above conditions, DCO + op- 
tically thin lines constrain the remaining nonthermal 
broadening (turbulence in the following). Satisfactory 
fits to the DCO+ J=3-2 line (with [DCO+/H 2 ] = (l- 
2)xl0 -10 ) are obtained for er tur ; ) =0.15kms~ 1 (the tur- 
bulence on the size scale of the beam). Therefore, 
(Jturb m the densest parts of the envelope (as traced 
by the DCO + ) is subsonic, and probably a vestige of 
the low degree of turbulence found in prestellar cores 
(Bcrgin & Tafalla 2007). Hereafter we adopt a uniform 
cturb value of 0.15kms _1 in the models of the envelope. 

In order to reproduce the HCO + J =1-0 line pro- 
file we have added the contribution from the lower 
density (more turbulent) ambient cloud surrounding 
IRAS 22129+7000 and producing a line-of-sight extinc- 
tion of Ay =2-3 (from 1.2 mm observations). Given the 
high HCO + dipole moment and large abundance (pro- 
ducing subthermal excitation and large opacities), low- 
J line photons from the dense envelope will be greatly 
attenuated and scattered over large areas by this low 
density halo (Ccrnicharo & Guelin 1987). Note that the 
negligible DCO + and C 18 column densities in the ambi- 
ent cloud do not alter the DCO + and C 18 line emission 
from the envelope itself. The density and temperature in 
the ambient cloud can be constrained with the observed 
CO lines. They are reproduced with n(H2)~1000cm -3 
and Tfc ~15K, in agreement with estimations by Witt et 



al. (1987). The halo ambient cloud in the model is at rest 
(uinf^Okms" 1 ) but the broader observed linewidths re- 
quire (Tt llr f,=0.30kms _1 (supersonic), at least a factor 2 
larger than in the circumstellar envelope. Figure 5 shows 
the best fits and the density and temperaure gradients 
used in the model. In spite of the good agreement, the 
model slightly underestimates the observed HCO + and 
C 18 J=l-0 linewidths, likely due to a minor contribu- 
tion from the high velocity outflow (not modelled here). 

The resulting large DCO + /HCO + abundance ratio 
(~0.02) found in IRAS 22129+7000 is similar to that 
derived in cold and dense prestellar condensations near 
other PDRs (e.g., Pety et al. 2007). This is consis- 
tent for a species that shows enhanced abundance in the 
cold FUV-shielded gas (low ionization degree) affected 
by ongoing deuterium fractionation (where dissociative 
recombinations with electrons are its main destruction 
route). As a consequence, the DCO+/HCO+ abundance 
ratio is roughly inversely proportional to the electrons 
abundance and we can estimate its value. Taking the 
same gas phase chemical network used by us in Pety 
et al. (2007), and assuming a cosmic ionization rate of 
C=3xl0 _17 s _1 (Dalgarno 2006), we derive an ionization 
fraction of ~10~ 8 (with respect to H nuclei) in the dens- 
est and coldest regions of the envelope where the DCO + 
abundance peaks. This value is, within a factor ~2, sim- 
ilar to the ionization fraction inferred in prestellar cores 
such as Barnard 68 (Maret & Bergin 2007). 

3.4.2. CO outflow parameters 

In the previous section we have estimated the amount 
of material towards IRAS 22129+7000 and also evalu- 
ated the CO line emission contribution from the ambient 
cloud and from the circumstellar envelope (through the 
optically thin C 18 line). After subtracting the CO con- 
tribution from the above line-of-sight components we 
now estimate a CO column density of ~10 16 cm -2 in the 
outflow. Taking into account the observed dimensions 
of the blue- and red-shifted lobes we compute the mass 
(~0.14M©), momentum (~1.OM kms" 1 ) and kinetic 
energy (~6.8xl0 43 erg) carried out by the outflow. Note 
that the mechanical luminosity of the flow (^0.02 L Q ) is 
still a significant fraction (~0.4%) of the radiant luminos- 
ity of the driving source. The estimated outflow momen- 
tum flux is F C o ~3xlO~ 5 M kms- 1 y r ~\ although the 
correction factor for inclination effects (sin i / cos 2 i) will 
decrease the actual Fco if i <40°. Hence, Fco is inter- 
mediate between the observed momentum fluxes found 
in early Class sources and those found in late Class I 
sources (e.g, Bontemps et al. 1996, Bachiller 1996). 
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TABLE 3 

IRAS 22129+700: CO outflow parameters 









Kinetic 


Mechanical 


Outflow 


Mass 


Momentum 


Energy 


Luminosity 


Component™ 


(M ) 


(Mq kms -1 ) 


(10 43 erg) 


(lO- 3 L Q ) 


Blucshiftcd... 


0.06 


0.4 


2.9 


6.9 


Redshifted.... 


0.08 


0.6 


3.9 


9.3 


Total 


0.14 


1.0 


6.8 


16.2 



a To calculate the energy and momentum we use an expansion velocity 
of 7kms _1 and a dynamical age of 35,000 yr (see text). 



4. DISCUSSION AND CONCLUSIONS 
4.1. The Nature of IRAS 22129+ 7000 

As noted earlier IRAS 22129+7000 shows indirect ev- 
idence of a protostar/disk system (a collimated out- 
flow) and it shows extended dust emission from the 
envelope with a large L smm /Lb ; ratio (~2%). These 
characteristics are consistent with the standard observa- 
tional definition of Class sources (Andre et al. 2000). 
Comparison of the derived Lf, Q ; and M™„ values with 
theoretical evolutionary tracks (Saraceno et al. 1996; 
Andre et al. 2000) shows that the source is likely in the 
interesting Class to I transition. Approximately 50% 
of the initial mass envelope must have been accreted at 
this stage. According to these diagrams, infall must have 
started in the parental core ~(0.5-l) xlO 5 yr ago. This 
age is consistent with the dynamical age of the outflow, 
specially if inclination effects are taken into account. 

The peculiar SED of IRAS 22129+7000 shows however 
a near- and mid-IR excess that produces a high Tb i 
of 175 K, consistent with more evolved Class I sources 
(~70-650K; Myers et al. 1998). Given the strong de- 
pendence of the short-wavelength IR emission on the 
viewing angle (Fig. 3), the inferred T{, ; is certainly af- 
fected by orientation effects. Up-to-date 2D SED models 
show however that inclination is not enough to repro- 
duce the observed excess. Nevertheless, as mentioned 
by Robitaille et al. (2006), their large grid of modelled 
SEDs have some caveats. The most important for our 
analysis is that none of them includes multiple central 
sources. An interesting possibility is thus that the source 
in IRAS 22129+7000 is a multiple (e.g., binary) stellar 
system. In fact, prestellar cores can fragment a second 
time during protostcllar collapse and lead to the forma- 
tion of binary systems (e.g. Duchenc ct al. 2007). How- 
ever, it is not yet clear if the initial conditions in these 
cores (e.g., the cloud environment) play a major role in 
such dynamical fragmentation. Simulations predict that 
external compression waves arriving at the central object 
at different times or instabilities in the outer accretion 
disk could be the seed of multiple protostellar systems 
(e.g., Hennebelle et al. 2003, 2004). In this picture, the 
presence of both a young and a more evolved YSO shar- 
ing a common circumstellar envelope would explain the 
observed near- and mid-IR excess if the scattered light 
from the evolved YSO escapes from the outflow cavities. 



Higher angular resolution IR observations are required to 
confirm this scenario. In addition, a larger set of molec- 
ular line observations and maps are needed to tightly 
constrain the physics (infall and rotation velocity fields) 
and chemistry (molecular content) of this protostar and 
its environment. 

4.2. Star Formation Near Ced 201 PDR 

IRAS 22129+7000 is located near Ced 201 PDR 
(<i^0.4pc and Ay^l). Taking into account geometri- 
cal dilution and selective dust extinction of FUV disso- 
ciating radiation from BD+69°1231, Casey (1991) com- 
puted that a region with a radius of at least ^0.6 pc 
(^300") around the exciting star is permeated by FUV 
photons. In particular the FUV radiation field near 
IRAS 22129+7000 should be presently at least -2 times 
the mean interstellar radiation field. Aside from dynam- 
ical effects associated with the presence of a nearby PDR 
illuminated by a runaway star (e.g., shock waves), heat- 
ing and radiation pressure in the region are larger than 
in other FUV radiation-free environments. Interestingly 
enough, the 1.2 mm dust emission (very weak in the PDR 
itself) peaks again in two dust condensations ^0.25 pc 
North of Ced 201, roughly in the same direction of out- 
flow red-lobe axis. These cores/clumps have a radius 
of ^5,000 AU and coincide with an enhancement of the 
HCO + J=l-0 line emission (see Fig.l). It is very likely 
that these condensations are being externally heated and 
compressed by the combined action of the outflow and 
the PDR associated FUV-driven shock wave. In fact, 
a similar scenario has presumably led to the formation 
of IRAS 22129+7000 itself. We thus suspect that star 
formation in the B175 globule is triggered by the shocks 
and radiative feedback induced by BD+69°1231 star and 
the new generation of low-mass YSOs. 

Observations of cores and YSOs provide new insights 
into how the star formation process is affected by its 
cloud environment. The vicinity of late type B stars has 
interesting intermediate properties between isolated dark 
cores (where FUV-radiation does not play a major role) 
and cores compressed by the expansion of H n regions 
(near more massive stars). Cool PDRs such as Ced 201 
are nearly devoid of ionizing photons and thus they pro- 
vide a very different environment compared to other star 
forming regions which are photoionized by brighter OB 
stars. In the special case of Ced 201, the exciting star 
is moving fast through the region, suggesting that its 
arrival to the cloud, and thus the triggering, is recent. 
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